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FOURTH  ORDER  EVOLUTION  EQUATIONS  FOR 
SHOALING  NAVE  GROUPS 

Ruth  Iusim  and  Michael  Stiassnie 
Department  of  Civil  Engineering,  Technion 
Haifa  32000,  Israel 

Third  order  evolution  equation  for  groups  of  water-waves  moving 
over  an  uneven  bottom  were  first  derived  by  Dj ordjevic'  and  Redekopp 
(1978).  In  the  present  report  we  carry  their  derivation  one  step 
further,  to  fourth  order,  utilizing  the  REDUCE  2  algebraic  manipulator, 
see-Heam  (1973). 

We  consider  the  evolution  of  a  uni -directional  progressive  gravity 
wave  moving  along  the  x  axis  on  the  free  surface  of  a  homogeneous  liquid 
with  depth  h  ■  h(x)  varying  in  the  direction  of  the  pr&pagation. 

The  effect  of  surface  tension  is  neglected,  so  the  analysis  applies  to 
gravity  waves  only.  The  fluid  motion  is  inrotational ,  thus  there  exists 
a  velocity  potential  $(x, z,t)  satisfying  Laplace's  equation 

0  <» 

where  z  is  the  vertical  coordinate,  and  z  ■  0  is  the  undisturbed 
free  surface. 

The  boundary  condition  on  the  bottom  z  *  -h(x)  is 

♦i  ■  -h'<x>  \ 
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where  9  *  /  k(e)dx  -wt 
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and  c.c  means  complex  conjugate. 


With  e  chosen  to  be  small,  the  functions  ^(£,t,  z)  and 
Cj(C,T),j  >  1  are  expanded  formally  in  powers  of  e  as  follows: 


*.(e,T,  z)  *  Z  e  *  -CC.t.z) 
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We  also  assume  ^(C.TjZ)  <  0(e) 
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thus  only  wave  induced  currents  are  considered. 

Substituting  the  Fourier  series  (6)  into  the  Laplace  equation 
(1)  we  obtain  that  the  zero  order  potential  satisfies: 

V2*0U,t,z)  *  0  (12) 

From  (11)  and  (12)  it  follows  that: 

♦c  (C,x,z)  *  0(e3)  (13) 

zz 

Expanding  the  free  surface  conditions  (3)  and  (4)  around  the 
equilibrium  position  z  ■  0,  substituting  eqs.  (6)  to  (10)  and  looking 
for  the  coefficient  of  e°,  yields  respectively: 
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and 
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where  *  ♦d(C,t,o),  ♦.j  *  ♦^j(C»T,o),  and  the  bar  denotes  the 
complex  conjugate. 

From  (14)  and  (15)  it  follows  that  the  order  of  CQ  is  greater 

or  equal  to  e2  and  that  the  order  of  4  (0)  is  greater  or  equal 

°z 

to  e3 


CQ  <  0(e2)  -  (16) 

and 

♦Q  (P)  «  0(c3)  (17) 
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The  next  step  of  the  derivation  is  similar  to  that  given  by 
Djordjevic*  and  Redekopp  with  the  difference  that  we  continue  the 
process  to  fourth  order.  Substitution  of  (6)  (8)  and  (9)  into 
Laplace's  equation  (1)  gives: 
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Substitution  of  eq  (6)  to  (10)  and  eqs.  (18)  to  (23)  into  the 
expansion  around  z  =  0  of  the  free  surface  conditions  (3)  and  (4) 
gives : 
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where  y1#...,y5  are  given  by  (32)  to  (36)  and 
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Combining  eqs.  (31)  and  (40)  and  introducing  the  new  variable 

Aj  *  eA  +  e2D  (5C 

one  finally  obtains  the  modification  to  fourth  order. of  the  cubic 
SchrtJdinger  equation 
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On  the  other  hand  substitution  of  (18)  and  (24)  into  (14)  and  (15) 
and  elimination  of  yields  the  following  equation  for  4>q: 
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which  can  be  written  to  the  same  order  of  accuracy 
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Laplace's  equation  (12)  and  the  bounadry  conditions  (2)  at  z*  -h 

and  eqs.  (51),  (53)  at  z  =  0  form  the  system  of  equations  from  which 

A.  and  4*  can  be  determined. 

1  o 

For  the  case  of  infinitely  deep  water,  the  terms  X  *e2  A.  <J>  , 

0  XT 

6  ,  and  X-e^.*  in  (53)  are  0(e5),  and  equations  (53)  and 

/  o  1 

(51)  become  respectively: 
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The  last  two  equations  are  identical  to  those  given  in  K.B.  Dysthe 
(1980). 


Note  that  the  following  typographical  errors  were  found  in  Dysthe's 

paper: 

His  equations (2. 17)  and  (2.19)  should  be  written  as  follows: 

r  =  4k4|Al2+  8ik3(AA*  -A*AJ“4ik3[A|2  +2u»k(f  -if  ) 

XX  X  X  z 

(Dysthe  2.17) 

2i<V’A)*!VV|A|2=  * 

+  |  iA(AA*  -A*Ax)  -h\  A|  2Ax  +  A($x-if2)  (Dysthe  2.19) 
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